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o : 

O ■ Abstract 

(N ; 

| The radial motion of matter in a centrally symmetric gravitational field in 

a comoving reference frame is investigated for a realistic equation of state of 
matter. The dynamics of the formation of an event horizon are investigated. 

r-> : 

^ : 1 INTRODUCTION 

<n : 

The formation of a black-hole event horizon has attracted a great deal of attention 
on the part of physicists for a long time. An enormous amount of material has 
been written on this subject (see, for example, Ref. [p], ||, |[ £|); nevertheless, the 
treatment of this problem within the general theory of relativity has created more 

*q \ questions than known solutions. 

One of the main questions concerning this problem is still the reciprocal influence 
of accreting matter on a black hole. The motion of test particles in the field of a 
black hole has been considered hitherto for the most part, but they, as we know, 
do not exhibit a reciprocal influence, which can be enormous when a falling particle 

^ | achieves the speed of light as it crosses the event horizon. 

1. In this paper we consider the special, but physically real case of spherically 
symmetric accretion on a central body without allowance for rotation. The following 
notation is adopted: the speed of light c and the gravitational constant G are set 
equal to unity. In these units the gravitational radius for a given mass M is r g = 2M, 
i.e., the radius of the event horizon in free space for the same mass concentrated at 
the center. 

Let us devise a likely model for the evolution of the system. We assume that 
our system is a cooling massive star having a radius Ro and a gravitational radius 
Rg = 2M(Rq), where Rg < Ro- The matter comprising this body is initially 
at rest ("dust" with the equation of state P = as, where P is the pressure in the 
matter, e is the energy density, and a is a constant). In the next moment the 
matter begins to fall freely. [|TJ ) If it is assumed that the gravitational fields are not 
excessively strong and that the dust density pj ) is fairly small in the initial moment, 
a force field with a finite energy is needed to retain it in the initial moment. After 
this field is removed, the dust leaves the system and ceases to interact with it after 
a time of the order of the size of the system, i.e., after a time much shorter than the 
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time during which the dust manages to partially settle and the gravitational fields 
increase dramatically. Thus, this model is physically consistent. 

What subsequently happens to the system? The dust begins to fall toward 
the center of the body, increasing its mean density and the gravitational radius 
r g (r) for the mass M(r) at a certain radius r. If we would neglect the reciprocal 
influence of the pressure of the moving matter on the dynamics of the system and 
on its gravitational field, then after all of the matter has unavoidably fallen and the 
inequality 

r g (r) = 2M(r) > r (1) 

holds at one of the points r of the system, an event horizon would form at that point 
according to Schwarzchild's solution for a gravitational field in a vacuum, i.e., the 
velocity of the falling matter relative to the r = const surfaces would reach the speed 
of light (see below). Is this what actually happens? The achievement of the speed 
of light by the matter causes a change in the sign of the interval and is therefore an 
invariant event, which does not depend on the choice of the reference frame. 

An attempt to solve this problem in a reference frame which is stationary at in- 
finity leads at once to a contradiction. In fact, the analytically exact, nonstationary 
model, in principle, cannot be studied. If, however, a simplification is made and it 
is assumed that the system is quasistationary at a certain moment in time in the 
range from one radius to a certain radius known to be large, but still far smaller 
than the dimensions of the system, then it can be stated, at the very least, that 
the components gttij) and g rr (t) of the metric have singularities (zeros and poles) 
in this reference frame. [ft]|p When the parameters of the system are chosen so that 
there would be a region in space where the inequality (1) is sure to be satisfied, it 
becomes clear that the metric does not have singularities, regardless of whether the 
inequality (1) is satisfied. 

This can be shown by assuming that if a singularity appears at a certain point 
To, the component of the metric near it can be represented in the form 

g u (r) const(r - r ) y \ 

where is a certain number. When such a metric is substituted into the equations, 
it is found that they do not have a solution for any yt ^ 0. 

This apparently indicates that the singularities and thus the horizon of the 
rapidly moving matter are eliminated (the right-hand side of the Einstein equations, 
which is equal to zero in a vacuum, becomes singular in the presence of ultrarela- 
tivistic falling matter when the radial component of the three- velocity tends to unity 
and the radial and temporal components of the four- velocity tend to infinity; this is 
also the reason for the elimination of the singularities of the metric). 

However, in reality all this stems from the inapplicability of the quasistationary 
approximation in the case of strong gravitational fields. It is inapplicable because 
the passage of time in the system is highly nonuniform due to the nonuniformity 
of the component <7#(r) of the metric. This causes the picture, which appears to 
be stationary far from the center, to become highly nonstationary to an observer 
approaching the symmetry center of the system. 

Nevertheless, this does not remove the question posed: do an horizon and a black 
hole appear in the real nonstationary case? 
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2. An answer to the question posed can be found by selecting a co moving 
reference frame. The problem was solved in this frame in Ref. 1 (Sec. 103) in the 
special case of a = (see below). Matter is at rest in the reference frame chosen, 
and its motion can be evaluated only from the variation of the "circumferential" 
or photometric distances r, which are related to the center of the system and are 
defined as the circumferences of the respective circles around the center: 27rr. When 
the radius r is defined as such, it is convenient to represent the metric in the form 

ds 2 = e u dt 2 - e x dR 2 - r 2 (d6 2 + sin 2 9 dip 2 ). (2) 

Here R is the coordinate of a dust particle in the comoving reference frame or its 
index, and e v , e A , and r are functions of R and t. It is noteworthy that at zero 
pressure, i.e., when a = 0, we have v — 0, i.e., the reference frame is simultaneously 
synchronous. 

To solve the problem posed we write out the Einstein equation in the comoving 
reference frame: 

r /2 e~ A (l + ru'/r') - e~ v {2rr + r 2 - rrv) = 1 + 8nar 2 e, (3a) 

2/1' + fifi' - V - u'fi = 0, (3b) 
2 \ 

X + 2/2+ Ins =0, (3c) 

1 + a J 

In \ ' 

v + — — lne) =0. (3d) 



1 + a 

Here jj, — 21nr, a prime denotes differentiation with respect to R, and a dot denotes 
differentiation with respect to t. Equations (3) were derived in Ref. [l| [Eqs. (2), 
(5), and (6) of problem 5 in Sec. 100]. 
It follows from (3d) that 

In 

~ \ne + f*(t) 



1 + n 



and that by transforming the time t in the interval element (2) the function f*(t) 
can be set equal to [2a/ (1 + a)} lne*, where is a constant with the dimensions of 
energy density, which expresses the measurement scale of e. Then 

v = In — . (4) 

1 + a £* 

We next assign the indices R to the dust particles so that r = R in the initial 
moment. Under such initial conditions r'(R,t) corresponds to (rio/n) 1 ^ 3 , where 
n(R, t) is the concentration of dust particles and no is its value at the initial moment. 

Let us now ascertain the conditions which must be imposed on the initial distri- 
bution of the dust. The most important among them is that the inequality (1) need 
not hold within the matter at the initial moment. It means that there is no horizon 
in all space in the initial moment. It thus imposes an upper limit on the initial 
density of the dust and on the initial dimensions of the system. More specifically, if 
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the initial density distribution of the dust is set equal to So(R), then, according to 
(1), the maximum radius of the body -R max is uniquely specified by the expression 



Rn 



Rr, 



Ane (R)R 2 dR. 



(5) 



Then, it follows from (3c) and (4) that 



d 



-KA + 2//)-H 



or 



u = a[X + 2fx + r(R)}, (6) 

where f*(R) is an arbitrary function that depends on the initial conditions. 

3. Let us now find the initial values for all the variables in our problem. We 
have already assigned these values for r and e. From (4) it follows that 



2a e 
In — . 

1 + a 



(7) 



To find the initial value of A we take advantage of the fact that the problem has 
already been solved for a = 0, and we can therefore utilize the familiar expression 
for AqU=o from Ref. (Sec. 103.6): 



\ (R) = -ln[l-S(R)}, 



(8) 
(9) 



where for a = we have 

S(R) = 2M(R)/R, 

and M(R) is the mass within the radius R at the initial moment. 

The expression for S(R) for an arbitrary value of a is the same. It can be 
obtained from Eq. (4) in problem 5 of Sec. 100 in Ref. |l[], where the Einstein 
equations in matter in the comoving reference frame were found for a centrally 
symmetric system. We write out this equation: 



— e 



—X 



2e v 



V+(A) 2 /2 



8tt£. 



(10) 



Expressing \i in terms of r (fi = lnr 2 ) and combining similar terms, we can bring 
this expression into the form 



Inr'sr 2 = —\r(r' 2 e A 



/ r * 

1)]' H Arf + (f 



Taking into account the expression (8), as well as the fact that, according to the 
expression (100.23) in Ref. M, the equality 



2M(r) 



lire(f,t)f 2 dr\ t=const 
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holds for the initial moment in time, when f = and r' — 1, we obtain the expression 
(9) for S(R) after preliminarily integrating (11) over R from to R. 

Substituting the expression (8) into (6) and taking into account (7), we find that 



f*(R) = —\ n ^ + \ n [i-S(R)]-\nR\ 

1 + a 



;i2) 



4. Now, plugging (6) into (3b) and dividing everything by fi/i', we obtain the 
expression 

-(21nu + u — A) — — — a . (13) 

Taking into account that e~ u (2rf + r 2 — rrv) = (e~ u rf 2 )' /f and introducing the 
notation 

U(R,t) = (r) 2 , Q(R,T) = r' 2 e-\ 
we see that Eq. (3a) can be written as an equation for U : 



(14) 



U 

h at/ = <7, 

r 



(15) 



where 



a(R,t) 



TV 

r 



a(R,t) 



rv 



Q 1 + — - 1 - 87rar 2 £ 



This equation has a solution which satisfies the initial conditions: 

t rt 
U(R,t)^ : ^—^l^(R,t)a{R,t)fdt, 1 *(R,t) = ex V J a(R,t) 



(16) 



Finding £/, we can obtain an expression for the square of the velocity of the matter 
relative to the r = const surfaces from the form of the metric (2) (see Appendix 1): 



(17) 



V 2 (R,t) = Ue- u e x /r' 2 . 

The expression for 7* can easily be found: 

j*(R,t) =C{R)re-\ 

where the multiplier C(R) for 7*(i2, t), which does not depend on t, can be taken 
out of the integral sign in (16) and canceled; therefore, it can be set equal to unity. 
Then 

1 + — J - 1 - Mar 2 e. 

Alternatively, taking into account that Eq. (13) can now be rewritten as 

(2 In // + /!- A)' (lnQ)' v' 



we obtain 



, (r(Q-l))- 2 
7 cr = : 8%ar s. 



(18) 
(19) 



Then (16) is rewritten in the form 



U = — [r(Q-l) -R(Q -1) + 2am(R,t)\, (20) 



where we have introduced the notation 

o R 



m(R,t) = J AiT£r 2 rdt = J 4ner 2 r dr\ R=cosnt . (21) 



t r 



5. For a = 0, taking into account (6), (9), (18), and (20), we can easily obtain 
an analytically exact expression for U and V: 



U a=0 = S(R) ( - - l) = - S(R). (22) 

V r / r 

Substituting this expression into (17), for the velocity we obtain 

2 2M(R)/r - S(R) 
K=o- 1 _ S (R) ■ W 

Hence V a= o = 1 when r = ro = 2M(R). This coincides with the results in Sec. 100 
of Ref. [l[], where the problem has already been solved for this case. 

6. Finally, let us consider the location of the horizon|TJ]p in the presence of a 
nonzero pressure. For this purpose we plug the expressions (20) and (8) for e~ x ° into 
formula (17). After some relatively simple transformations, we ultimately obtain 

2M(R)+2am(R,t) 
r = ' < 24 > 

As will be shown in Appendix 1, the horizon appears at the point and at the time 
where the velocity of the falling matter relative to the r = const surfaces reaches 
unity, i.e., where V — 1. In addition, the speed of light relative to the falling matter 
at this site is also, as always, equal to unity. 

Hence, according to (24) and (9), the horizon radius r^ or is given by the formula 

r hor = 2M(R) + 2am(R, t hoT ) . (25) 

Thus, the horizon is displaced to a larger radius in comparison to the value in a 
vacuum r 0hor = 2M(R) by 2am(i?, t hor ). In this case the quantity m(R,t) has the 
meaning of the mass which would accumulate if we would join layers of dust with 
the initial radius R and the thickness dr(t) to one another up to the radius r(R,t) 
at the moment when this dr layer passes through the joining point. 

7. Regarding the possible values of a we note that a = corresponds to dustlike 
matter without interactions between the particles. The results obtained for them 
are the same [see (23)] as the results for test particles in a central field of mass 
M (see Sec. 101 in Ref. |]|). However, of course, such an equation of state of 
matter cannot correspond to reality near the horizon. It is reasonable to assume 
that the ultrarelativistic equation of state of matter, in which a = 1/3, holds near 
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the horizon. Therefore, the location of the horizon should probably be sought with 
just such a value of a. 

8. When a ^ 0, it would appear that the falling matter should be slowed under 
the action of the pressure gradient, and the horizon should therefore form later, i.e., 
be displaced toward smaller values of r, but, as we have just shown, it is displaced 
toward larger values of r by 2am(R, th or ). What is the reason for this contradiction? 
It can be seen from the initial equations (3) that the reason should be sought in Eq. 
(3a). For this purpose we explore Eqs. (3a) and (4) in the initial moment for the 
case of a <C 1. In that moment r = and r' = 1; therefore, we write 

[1 - S(R)\ (l - 2ar^ - 2rf {^j^ ~ 1 + 8nar 2 e. 



Since 



(e/e*) 2a *l + 2a\n-, S(R) = r 
£* R 



R, 



then, after performing some relatively simple transformations, in the linear approx- 
imation with respect to a we obtain 



GM(R) 





VP 


1 - 2a In — 




L £*. 


P 



1 - 



2GM(r) 



rc z 



— AnaGrp, 



(26) 



where p(R,r) = e(R,r)/c 2 is the density of the matter. Here, for the sake of clarity 
we use the ordinary (Gaussian) system of units with G ^ 1 and c ^ 1. It can 
be seen from (26) that the first term corresponds to the ordinary Newtonian force 
of gravity, and the second term corresponds to the interaction force between the 
particles, i.e, the pressure gradient (just this force is the cause of the slowing of 
the fall of the matter in the first stage). The remaining terms do not appear in 
the equation of motion in the Newtonian approximation (the corrections in square 
brackets are also neglected in that case), but, as we have already seen, the last 
term begins to dominate over the second term at high energies; therefore, a shift of 
the horizon toward larger radii appears. Thus, the contradiction has been resolved. 
Physically this corresponds to the "gravity of pressure" in the general theory of 
relativity, which surpasses the gradient terms at high energies. 

9. The analysis performed allows us to draw the following conclusions. 

First, a shift of the horizon toward a larger radius in comparison to the Schwarzchild 
radius due to the "gravity of pressure" has been discovered. We stress that this ef- 
fect is purely dynamic and is not observed in the static case (after all the matter 
has fallen). 

Second, according to the results in Appendix 2, the evolution of the entire system 
at a constant value of a is completely specified by the energy density distribution 
profile in the initial moment, i.e., for example, by the normalized density distribution 
of the matter and by the value of the parameter S at an arbitrary point on this 
distribution. 

If the evolution of only one spherical layer of matter with the index R must be 
described, it is completely specified by three dimensionless parameters in the initial 
moment in that layer and, in this sense, does not depend on the initial distribution 
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of the matter in the system below and above that layer. However, this in no way 
signifies the independence of the spherical layers in the general case, since just these 
three parameters, as will be seen from Appendix 2, govern the interaction of the 
layers. Consequently, integration of the system leads to a complete family of self- 
similar solutions. 

Third, according to Appendix 2, a local extremum appears on the V(R)\ t = C onst 
curve for a specific choice of initial parameters, and when V — 1, it leads to the 
formation of a second apparent horizon in the system (an analog of the second 
horizon in the Reissner-Nordstrom and Kerr-Newman solutions for an electrically 
charged rotating static black hole; for an interpretation of these solutions, see, for 
example, Refs. || and ||). 



We have hitherto used the term horizon to refer to a trapping surface, or an apparent 
horizon, as it is called in the literature. 

Let us ascertain the difference between an event horizon and an apparent horizon 
in greater detail in an example. We assume that we already have a stationary black 
hole of mass M and that there is an apparent horizon at r = 2M. Now we assume 
that another chunk of matter with a mass 5M falls into our black hole. After it falls, 
the radius of the apparent horizon increases to 2(M + 8M). Thus, if an observer 
is placed between these radii before the additional chunk of matter falls, he would 
then be outside the black hole, but after the chunk of matter falls he would be inside 
it. The concept of an event horizon is global and is determined by the entire course 
of evolution of the black hole or, stated differently, by all the mass which falls into 
it at any time. 

The existence of an apparent horizon, which specifies a black hole locally, is 
sufficient for the existence of a black hole. As follows from our arguments, in the 
spherically symmetric case the two horizons ultimately coincide and form a static 
black hole described by Schwarzchild's solution. Therefore, we shall henceforth use 
the term horizon to refer to the apparent horizon. 

Let us prove that the horizon in a system with spherical symmetry forms at the 
moment when a falling particle with a nonzero rest mass achieves the speed of light 
relative to the r = const surfaces at the same point. For this purpose we write the 
law of motion for the particle in the form 



We now assume that we are located on a dust particle with the index and we 
are tracking a dust particle with the index R p , which sends us a light beam passing 
through the radii r p (R p , t), from the large radius r^t). The criterion for determining 
that the dust particle has not yet reached the horizon is the fact that we still see 
light from it, i.e., the light propagating still crosses the radii r > r p . Therefore, the 
criterion for determining that the dust particle has reached the horizon is an event 



1 




(27) 
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Figure 1: 

in which the light propagating from R p can no longer cross the radii r > r p . Let us 
express this criterion mathematically. 

In Fig. [l| the vertical straight line abode denotes the world line of an Rp dust 
particle in the coordinates R and t of the comoving reference frame from the moment 
of rest (a) to the center of the system (e) at r = 0. In this case of solid curves passing 
through points e, d, c, and b denote, respectively, lines of constant values of r(R,t) 
for r = 0, r < r hor , r = r hor , and r > r hor . The dashed lines emerging from these 
points denote the cones within which light emitted by the R p dust particle can 
propagate (light cannot propagate outside these cones). Therefore, according to the 
criterion indicated above, the horizon forms at the point where the cone is tangent 
to the r = const line. In the figure this line is designated as r = r hor , and it passes 
through point c. For clarity, Fig. |] shows that the light cone intersects lines with 
r > r p at point b; therefore, there is still no horizon at that point. This figure also 
shows that at point d the light cone is located entirely above the r = const curve 
passing through point d. Consequently, this light cone intersects only lines with 
r < r p , and therefore point d is already located below the horizon. 

Let us examine the expression (27) on one of the r = const curves and take its 
complete differential on that curve: 

1 f U'(R t) 
= dR-VUdt-- . 1 ' J dRdt, 

2 { y/U(R,t) 
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or 



Vu 



dt 

Jr 



i - 



U'(R,t) 
U(R,t) 



dt. 



(28) 



Next, differentiating (27) with respect to R, we obtain the following expression for 



r': 



t 



U' 



Vu 



dt, 



with consideration of which from (28) we find 

dt 
dR 



r=const 



(29) 



Thus, we have found an expression for the slope of an r = const curve relative to 
the R axis. 

To find the slope of a light cone, by definition, for light we have ds 2 = 0. Hence, 
from (2) it follows that 

' (30) 



dR 



Ve x - U . 



ligth 



According to the foregoing statements, the criterion for the absence of a horizon is 
the condition 

dt 
< dR 



dt 
dR 



ngth 



(31) 



r=const 



Substituting the expressions (29) and (30) therein and taking into account (17), we 
obtain this criterion in the form 



W\ 



< 1. 



(32) 



Here, according to (29), the rate of motion of the matter relative to the r = const 
lines has the form 



W\ 



dl_ 

dr 



= Ve x ~ u 



dt 
dR 



r=const 

Thus, the assertion that a horizon forms at the moment when the matter achieves 
the velocity V = 1 relative to the r = const surfaces has been proved. The horizon 
surface separates regions in which r is space-similar and time-similar. 



To solve the equations describing collapse, we first bring them into dimensionless 
forms. For this purpose it is convenient to introduce the following notation: 

_ p (R) _ 8ne (R)R 2 
' ' 7 (P) M(R) ■ 
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In this Appendix we find the ranges of permissible values of 7 and S, investigate 
the character of collapse at these values of the parameters, and obtain numerical 
solutions for V 2 . We must first of all know the form of the function r'(x). Differen- 
tiating (27), we obtain [pip 



1 x 

r\ x ) = 1 + -RJ [In U{R, x))' dx. 



(33) 



Unfortunately, an analytically exact expression for r' can be found only in the case 
of a = 0, the character of collapse can be assessed exactly only at that value of 
a. However, the main features of that character, as will be seen below, remain the 
same as in the case of a 7^ 0. Therefore, let us first investigate the case of a = 0. 
Thus, we should find r'(R, x). According to the expression (22) for U, we obtain 

\n[U(R,x)] = \n[S(R)} + ln f~ - lV 

Introducing the notation y = r' — x and taking into account that x' = y/R, we have 

The substitution of this expression into (33) gives 



y(x) + x- 1 = - 



RS' 
S 



y[x) 



x(l — x] 



dx. 



(34) 



Differentiating (34) with respect to x, we obtain 

dy(x) 



dx 



+ a*(x)y(x) = a*(R), 



(35) 



where we have introduced the notation 

1 



a*(x) 



2x(a — x) 



, tr*(R) 



RS' 
~2S 



— h -7. 

2 2 1 



As can be seen, Eq. (35) coincides in form with Eq. (15), and the initial conditions, 
y\ t= o = 0, are the same; therefore, the method used to solve it is similar. The 
solution has the form 



y{x) 



r — x = a 



1 — x 1 — x 
arctan \ / 



x 



x 



(1-x) 



(36) 



Let us find the domain of r'. First, the condition for compression of the matter 
has the form r' < 1. Second, the condition that dust layers with different R do not 
intersect[ 67 P has the form r' > 0. Thus, 



< r' < 1. 



(37) 
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We assume that the V 2 (R) curve for t — t m — const has a local extremum, and we 
presume (to fix ideas) that it is a maximum. Then the horizon appears specifically 
at the local maximum, i.e., the point R = -R ex tr- We now find the condition for a 
maximum. First, at that point we should have V 2 (R extr ,x) = V 2 Ktr . Second, since 
it is the first point at which the velocity of the matter achieves the value \4 xtr and 
the rate of collapse increases with time, in the vicinity of this point we should have 
r2 or 

extr' U1 

dV 2 (R,t m ) 



V 2 < V 2 



dR 



dV 2 (R,t r 
dR 



> 0, R < -Rextr, 



< 0, R> -Rextr- 



(38) 



If it turns out that (38) holds with opposite inequality signs, there will be a local 
minimum on the V 2 (R) curve at the point -R cx tr at the moment when the velocity 
Kxtr is achieved at that point, i.e., the matter will achieve the velocity \4xtr last at 
that point. 

The condition for an extremum is written in the form 



where, according to (23), 



dV 2 (R,x) 
dR 



V 2 (R,x)\ a=0 = 



0, 



1 - 1/x 



l-a/S(R)' 

Differentiating this expression with respect to R, we obtain 



dV< 



dR 



S/R 



t tr. 



-y + 



37 



(39) 



where it has been taken into account that x' — y/R and S'/S = (37 — 1)/R- Then, 
with allowance for the fact that < x < 1, < 7 < 1, —1 < y < 0, and < S < 1, 
the condition (38) can be rewritten in the form 



-y> 



-y < 



37-1 



37-1 



R< R, 



cxtr j 



(40) 



R> R 



cxtr • 



If we introduce the notation z — y(l — x)/x and take into account that, according 
to formula (23), z = Vy/l/S- 1, from (40) we obtain 



;(l-7) 



z arctan(^) 



1 + z 2 



37-I 
1-S 



> 0, R < -Rextr, 
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5(i 

2 V 



7) 



z arctan(z) 



z 



1 + z 2 

or for the extremum point we can write 

_,2 



37-1 
1-S 



< 0, R > R 



extr ; 



(1-7) 



zarctan(z) — 



1 + z 2 



S 



4/3 
1-S 



(41) 

This formula can be used to construct the plot of 7(5, V extr ) separating positive 
and negative values of the derivative {V 2 )' and to determine the character of the 
extremum. The corresponding curves for various values of V ext r are shown in Fig. 
[j]. The regions where V > are located above and to the right of them, and the 
regions where V < are located below and to the left of them. 

It is seen from Fig. [l| that there can be (for a definite choice of the distribution 
profile of the matter in the initial moment and of the parameter S at a certain point 
i?*) two values of R, at which V — 1 at a certain moment in time, and, therefore, 
the appearance of a second horizon in the system is possible. 

The appearance of a second horizon is not news in the physics of black holes 
(see, for example, the Reissner-Nordstrom or Kerr-Newman solution in Ref. J7[|). 

The results obtained in this Appendix apply to the case of the absence of pressure, 
although the case of a = 1/3 is of experimental interest. Therefore, we used formulas 
(4), (6), (8), (9), (12), (17), (18), (20), and (33) to introduce new dimensionless 
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variables (0 = u — u , X = X — X , U = Ue u °, and e = e/eq) and equations for them. 
The initial conditions for them take the form 

*>o = A = U Q = 0, r' Q = i = I. 

Designating the new coordinates as x = r/R and £ = R/R* (R* = const) and 



introducing the parameters fTJ 



h = £dtVo, h cv = S— 



1 + 3«7 



r] = h/h c 



we obtain equations for the new variables in the form 



e 0/a = x 4 e X 



-v/a\ (l+a)/2 



ln(r /2 e- A ) 



dtu h 

tytl tyt< 



dx, 



U = e v 



r' 2 e-\l-S)-l + --a^- [ ex 2 dx 
x x J 
1 



(42) 



r> = l + \ hd^U)di-hl-x). 



In the new variables the velocity is 



V 2 



Ue 



Hence Xhor = S + 2am, where 



i 

TJX f 

m = — = EqR 2 I An ex 2 dx, EqR 2 
R J 



3 7 g 



(43) 



We note that this formula and formula (36) can be used to find the corrections 
Srh OT in (25) to the displacement of the horizon in the linear approximation with 
respect to a, since, according to Eq. (103.11) from Ref. M for a = 0, we have 



8irEr z 



2M' 8tte R 2 



r' 



or 



ex 2 — l/r'. 



This expression can be substituted into (43) and a quadrature expression can be 
obtained for the correction sought. 
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Dependences of V 2 (x) at S=0.8. 
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Figure 3: 

In addition, we numerically integrated the equations for the case of a ^ using 
a difference scheme, and the results for various values of a are presented in Fig. [I]. 
As it should be, according to (25), the plots of V 2 (x) are displaced upward and to 
the right as a is increased from a = to a = 1/3. 

The numerical calculations confirm that the analytical results of this Appendix 
remain valid for the real equation of state of matter: P = as. Figure [I] shows 
plots of 7(5"), £(S), and rj(S) for the special case of a Gaussian density distribution: 
£ o(£)/ £ o(0) = exp(— 3£ 2 ). Comparing this figure with Fig. [l|, we can see that the 
7(5') curve in Fig. |1] crosses the 7(5) curves in Fig. [1| in the downward direction 
roughly at the point S ~ 0.92, if we proceed from £ = to £ = 1. As can be seen in 
Fig. [I], the point S ~ 0.92 corresponds to £ ~ 0.85 and rj ~ 0.1; therefore, since the 
region where V > is located above and to the right of the curves in Fig. [I] and 
the region where V < is located below and to the left of these curves, the point 
£ ~ 0.85 should be a local maximum on the V(£) curve for a constant value of t. 

This analytical result is confirmed by a numerical calculation of V(£)|i= CO nst 
curves, whose results are shown in Fig. ?? with the predicted maxima. 

To conclude this Appendix we would like to say a few words regarding the initial 
characteristics and distribution of the matter. 

When the equations of the model were brought into dimensionless form, it was 
found that the solution for a spherical layer of matter with the index R is com- 
pletely specified by three dimensionless parameters in the initial moment in that 
layer: 0<S'<1,0<7<1, and < r\ < 1. This corresponds to assigning the 
initial conditions for the gravitational potential and two parameters which deter- 
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Figure 4: 



mine the distribution of the matter and the pressure gradient near the point under 
consideration. Thus, upon integration we at once find a whole family of self-similar 
solutions, |8j]p which can be characterized by these three parameters alone and which 
contains the dependence on the other layers of matter above and below the radius 
R considered. 
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